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Abstract 

This  work  presents  new  variants  of  FFT  algorithms  (not  necessarily  power  of  two 
lengths),  in  a  computational  Kronecker-core  array  algebra  setting,  tailored  to  the  efficient 
point  spread  function  (PSF)  processing  in  multidimensional  synthetic  aperture  radar 
(SAR)  systems.  The  tailoring  of  the  algorithms  is  performed  through  a  targeted  study  of 
group  theoretic  properties  of  input/output  data  indexing  sets  and  associated  groups  of 
stride  permutations.  The  advantage  of  these  new  FFT  variants  over  conventional 
formulations  is  that  additive  group  theoretic  properties  of  multidimensional  input/output 
indexing  sets  are  used  for  their  mathematical  fonnulations,  establishing  mapping 
identifications  between  computing  structures  and  mathematical  expressions  identified  as 
factored  compositions  of  functional  primitives,  reducing  in  this  manner  their 
computational  complexity  and  improving  their  implementation  perfonnance. 

There  exist  many  fonnulations  of  fast  algorithms  for  computing  the  multidimensional 
discrete  Fourier  transfonn  (DFT).  In  this  work,  computational  Kronecker-core  array 
(CKA)  algebra,  a  branch  of  finite  dimensional  multi-linear  algebra,  is  used  as  a  language 
to  identify  similarities  and  differences  among  various  FFT  algorithm  variants  as  well  as 
for  the  creation  of  new  variants.  Each  multidimensional  DFT  computation  is  expressed  in 
matrix  form.  The  multidimensional  DFT  matrix,  in  turn,  is  decomposed  into  a  set  of 
factors,  called  functional  primitives,  which  are  individually  identified  with  an  underlying 
software/hardware  computational  constructs.  It  is  in  this  identification  process  where  the 
language  of  CKA  algebra  becomes  instrumental.  For  a  given  hardware  computational 
structure  and  multidimensional  DFT  matrix,  there  are  many  FFT  algorithm  variants 
which  can  map  to  this  target  machine.  The  language  of  CKA  algebra  aids  the  application 
developer  in  the  mapping  effort  of  identifying  the  more  computational  efficient  FFT 
variants  and  thus  reducing  the  computational  effort. 
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The  FFT  computing  methods  presented  are  currently  being  placed  in  a  Weyl-FIeisenberg 
computational  framework  setting  for  the  study  of  real-time  adaptive  computation  of 
linear  arrays  of  cross-ambiguity  functions,  formulated  as  point  surface  response 
functions  for  coded  and/or  linearly  FM  modulated,  time-frequency  collocated,  synthetic 
aperture  radar  (SAR)  transmitted  signals  and  received  signals  for  applications  in  signal 
detection,  information  extraction,  and  clutters  and  jammers  suppression  operations.  The 
Figure  One  below  depicts  a  basic  FFT-core  functional  primitive  developed  in 
MATLAB’s  Simulink  for  multi-FPGA  implementation  of  multidimensional  FFTs. 


Figure  One:  Basic  Core-FFT  Simulation 
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